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Abstract. The energies of the most energetic exten- 
sive air showers observed at the Yakutsk array have 
been estimated with help of the all detectors readings 
instead of using of the standard procedure with a 
parameter s(600). These detector readings have been 
compared with the detector responses, calculated for 
all particles which hit the scintillation detectors in 
each individual shower with observed the zenith and 
azimuth angles with help of the GEANT4 code. In 
turn the code CORSICA-6.616 have been used to 
produce such particles in the atmosphere in these 
individual showers and propagate them to detectors 
at the level of observation. Calculations have been 
carried out in terms of the QGSJET-2 and Gheisha- 
2002 models with the thinning parameter 10^ for 
the primary protons and helium, oxygen and iron 
nuclei. The energy of the most energetic extensive 
air shower observed at the Yakutsk array happened 
to be 200, 200, 180 and 165 EeV with the values of 
the function per one degree of freedom 0.9, 1., 0.9 
and 1.1 for the primary protons and helium, oxygen 
and iron nuclei accordingly. Thus interpreting data 
in terms of the QGSJET-2 and Gheisha 2002 models 
we conclude that after the bump and decreasing 
of the flux of the primary particles due to the 
Greisen, Zatsepin and Kuzmin effect that has been 
observed at the HiRes array and at the Pierre Auger 
observatory there is a possible feature in the energy 
spectrum - some increase in the flux at energies 200 
- 300 EeV observed both at the Yakutsk array and 
at the AGASA array. Such possible feature may be 
understood as the flux of heavy primary nuclei. It is 
also not excluded that some new component of the 
spectrum is observed or the Lorentz invariance may 
be violated at such huge energies. As an alternative 
conclusion the models QGSJET-2 and Gheisha 2002 
should be changed so as to produce much more 
muons in a shower. 
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I. Introduction 

The energy of extensive air showers (EAS) observed 
at the Yakutsk array (YA) is usually estimated as follows. 
First, for each individual shower the zenith and azimuth 



angles and coordinates of axis are determined. Then the 
signal s(600) is estimated from all scintillation detec- 
tor readings. This signal is determined by the energy 
deposited in a detector at a distance of 600 m from 
the shower core by all shower particles (by electrons, 
positrons, gammas and muons). Then this signal is recal- 
culated to the value it would have in the vertical shower 
with the help of the average value of the attenuation 
length < A > which is estimated by applying the 
constant integral intensity cut method |[T| . The calculated 
zenith-angle dependence is also may be used but with 
average value of the attenuation length < A >. But 
such procedure leads to large uncertainties in energy 
estimate for the individual showers. First, in terms of 
model calculation it was shown that the values of the 
attenuation length < A > estimated with the help of 
such procedure and calculated for the average shower 
development differ considerably |2 |, [3]. Second, and it 
is a main source of uncertainty, the individual showers 
may be generated by varies species of the primary 
particles. Besides, the real value of the attenuation length 
A in individual vertical showers may differ from ~ 200 
up to 2000 g ■ cm^^ |4| due to fluctuations in the 
longitudinal development. Indeed, the steepness of the 
individual cascade curve at the level of observation (and 
the attenuation length A) may vary very considerably 
due to fluctuations in the points of first and subsequent 
interactions of the primary particle with atomic nuclei 
in the atmocphere. Thus, the uncertainty of the signal 
s(600) estimated for the vertical showers with the help 
of a factor such as exp(— Ax/A) may be very large, 
where Ax is a slant depth. At last, the energy E is 
estimated with help of the following formula 



E = 4.8- W • s(600), eV. 



(1) 



In the formula (1) uncertainties are missed for simplicity 
and it was suggested that the signal s(600) is propor- 
tional to the energy E of a shower This formula is 
based on the calibration of signals with help of the 
Vavilov-Cherenkov radiation of a shower Again, this 
calibration has been carried out for the average value 
of a signal from some sample of showers. For the 
individual vertical showers the numerical coefficient in 
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(1) should vary due to fluctuations. Besides, it should 
be mentioned that in terms of the model calculation 
the value of this numerical coefficient was estimated 
as being 1.6 - 1.7 times less than quoted in (1) |5|. 
Thus, any alternative methods of energy estimation are 
of interest. It was suggested that readings of all detectors 
should be compared with calculated signals for a shower 
with the given values of the zenith and azimuth angles 
161. Calculations have been carried out for the giant 
shower observed at the YA fTl on the base of the original 
code. It was assumed in accordance with experimental 
data that this shower consist mainly of muons and their 
deflections in the geomagnetic field have been taken into 
account. The energy of this shower has been estimated 
as 3 • 102*^ eV. In this paper calculations of signals for 
this giant shower at many points with different distances 
from the shower core have been carried out for some 
sample of the individual showers induced by various 
primary particles to take into account fluctuations in 
the longitudinal and lateral development. Then the 
method has been used to find out which of calculated 
individual showers agree best with data fTj. As new 
energy estimate are happened to be rather high some 
analyses of energy spectra observed at various arrays 
have been carried out and the new interpretation of the 
energy spectrum has been suggested at ultra high energy 
region with the possible variable contribution from the 
local sources. 

II. Method of simulations 

Simulations of the individual shower development in 
the atmosphere have been carried out with the help 
of code CORSIKA-6.616 H in terms of flie models 
QGSJET2 |9| and Gheisha 2002 fTOl with the weight 
parameter e = 10^® (thinning). The program GEANT4 
ifTTll has been used to estimate signals in the scintillation 
detectors from electrons, positrons, gammas and muons. 
The bank of detector responces has been calculated for 
electrons, positrons and gammas with energies in the 
interval 0.001-10 GeV and muons with energies in the 
interval 0.3-1000 GeV which hit a detector at various 
the zenith angles (from 0° up to 60°). This bank of 
detector responces was used to estimate a signal in the 
scintillation detector when a shower particle hits it. The 
total area of 5 x 5 km? in the detector plane was divided 
into 201 X 201 squares with the side of 25 m. With the 
help of the code CORSIKA-6.616 the spread of shower 
particles in the detector plane has been estimated and the 
bank of detector responces has been used to calculate 
the signals in each square, regarded as a detector Thus, 
the matrix of 201 x 201 detector responces for each 
individual shower has been calculated. These matrixes 
of detector responces were calculated for individual 
showers with the same energy 10^° eV. Calculations 
have been carried out for four species of the primary 
particles (protons and nuclei of helium, oxygen and 
iron) with a statistics of four individual events for every 
species of the primaries. Readings of the 31 scintillation 



detectors have been used to search for the minimum 
of the function in th^ square with the width of 
400 m and a center determined by data with a step of 
1 m. These readings have been compared with calculated 
responces which were multiplied by the coefficient C. 
This coefficient changed from 0.1 up to 4.5 with a step 
of 0. 1 . Thus, it was assumed, that the energy of a shower 
and signals in the scintillation detectors are proportional 
to each other in some small interval. New estimates of 
energy, coordinates of axis and values of the function 
have been obtained for each individual shower 
The analysis of the energy spectra observed at various 
arrays has been carried out in the following way. The 
base universal spectrum Jb{E) — A ■ E^^-^^ has been 
suggested mainly on assumption of data lfT2]| at energies 
above 10^^ eV with A w 7.1 • 10^^ m-^s-^sr-^eV^-^^ . 
All possible features of the energy spectrum of the 
primary particles are considered relatively to this spec- 
trum. Besides, the reference spectrum Jr{E) has been 
suggested as follows. For the energy we will use besides 
E (in eV) additional notation y = lg(£'/l eV). In four 
energy intervals (i=l, 2, 3 and 4) 17. < y < 18.65, 
18.65 <y < 19.75, 19.75 <y < 20.01 and y > 20.01 
the spectrum Jr{E) has been approximated by the 
following exponent functions 



Ji{E) = 


A- 


E- 


-3.25 


J2{E) - 


C- 


E- 


-2.81 


UE) = 


D 


■ E 


-5.1 



J4{E) = Ji(S) = ^-^-3.25 

accordingly. Constants C and D may be expressed 
through A and equations for Jr{E) at the boundary 
points. For these four intervals we assume the reference 
spectrum as 

\gz,^lgiME)/J,iE)), 

where i=l, 2, 3, 4. This reference spectrum is then 
represented as follows 

Ig^i -0, 
lgZ2 = 0.44 -(y- 18.65), 
Ig Z3 = 0.484 - 1.85 • (?/ - 19.75), 

lgZ4 = 

accordingly. We consider the spectrum Jb{E) = A ■ 
^-3.25 universal up to highest energies. The first 
feature of the spectrum is suggested to be considered 
as some excess at energies 18.65 < y < 20.01. The left 
side of this exess is approximated as J2{E) — Ji{E) 
while the right side as J-i{E) — Ji{E). The possible 
second feature at energies (2-3)-102" eV win also be 
discussed. Results of the spectra J{E) observed at 
various arrays have been expressed as 

\gz = lg{J{E)/ME)) 

and are shown in comparison with the reference spec- 
trum. 
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Fig. 1: The energy spectra observed at various arrays 
and the reference spectrum (soUd hne): a) - lfT2ll . (b) - 
ISl, (c) - ISl, (d) - El 



III. Results of energy estimations and 

ANALYSIS OF THE ENERGY SPECTRUM 

The 16 various values of energy estimates for 16 
individual simulated showers with different values of 
the function have been obtained for the same sample 
of the 31 experimental readings of the observed giant 
shower Q. As the giant shower is very inclined muons 
make main contribution to signals in the scintillation 
detectors. The energy estimates are minimal and change 
inside the interval (1.6-1.75)-10^'' eV with the value 
~ 1.1 of the function per one degree of freedom 
for the iron nuclei primaries which produce more muons 
than the proton primaries . For the proton and helium nu- 
clei primaries energy estimates are maximal and change 
inside the interval (1.8-2.4)10^" eV with the value 
~ 0.9 of the function per one degree of freedom. 
For the oxygen nuclei primaries the energy estimates 
are in the interval (1.8-2)-10^" eV which is between 



intervals for proton and iron nuclei primaries. It should 
be mentioned that only muons contribute ^ 80% of the 
total signal and therefore energy estimate increases up to 
2.8 • 10^° eV for the proton primaries if only muons are 
taken into account as in |7| disregarding the contribution 
of electrons, positrons and gammas. New coordinates of 
shower axis vary from the experimental one by some 
dozen of meters. So we can make main conclusion that 
in terms of the QGSJET2 ID and Gheisha 2002 ||TOl 
models and the proton primaries there are extensive air 
showers with energies ^ 2 • 10^° eV. This estimate 
may be decreased up to ^ 1 • 10^° eV only in terms 
of some new model which produces twice as much 
muons as [9| and 1 10]. It is not easy because muons 
are produced in low energy region where the physics of 
interactions is known. But the energy estimate decreases 
up to ^ 1.6 • 10^" eV for the iron nuclei as the primary 
particles. It should be remind that such giant showers 
have been also observed in ITSl . It looks as a some 
contradiction to the suggestion by Greisen, Zatsepin 
and Kuzmin (GZK) lUH, ||T5]. Just on the contrary, 
the observations ||T21 . ||T61 show no giant showers in 
accordance with the GZK prediction but not with data 
|[T3l . iflTl . It is evidently that all world data should 
be understood. Of course, some uncertainties in energy 
estimates may exist. 

But it is also worth-while, to consider some new 
idea about the energy spectrum at ultra high energies. 
Usually, this spectrum is considered as universal and 
stationary [18|. It is possible if many uniformly dis- 
tributed sources contribute to the spectrum. In case of 
the near-by sources |19| distributed anisotropic-ally their 
contribution to the different intervals of the spectrum 
may be variable. If a power of local sources is not high 
one source may give a contribution to the number of 
observed events only once in many years. So, various 
sources may contribute to the spectrum at different time 
of an exposure. Due to deflection in magnetic fields 
arrival directions of showers differ from directions to 
local sources. If the number of sources is not high their 
contribution to the spectrum would be very variable 
and even chaotic. The primary particles from such local 
sources may be considered as some variable excess 
above the suggested universal spectrum Jb{E). Decreas- 
ing of the flux of the primary particles due to the GZK 
effect should probably be considered relatively to this 
universal spectrum but not relatively the observed bump 
as in |[T2l . lfT6l . Evidently, it is a difficult problem. 
Of course, the flux of particles at exess energies is 
also supressed by the GZK effect. Fig. 1 illustrates our 
suggestion. Data \gz = \g{J{E)/ Ji{E)) observed at 
various arrays are shown in Fig. 1 as follows: (a) - lfT2l 
(open circles - HiRes2, solid circles - HiResl), (b) - 
1 16 1 (solid circles), (c) - lfT3l (solid triangles) and (d) 
- 1 17 1 (solid pentagons). The reference spectrum is also 
shown on all Figures. As it was expected the data lfT2l 
agree very well with the reference spectrum. The data 
II 161 are below this reference spectrum, and the data lfT3l 
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and IfTTl are much above it. What possible features are 
seen in Fig. 1? First, no predicted dip |18| at energy 
10^^ eV is seen. Just the contrary is true: in this energy 
region the exponents of all spectra are approximately 
the same with a good accuracy. Second, all data show 
some excess in the energy interval (5-10)-10^^ eV (so 
called "bump" fTSl) which may be considered as a 
contribution from the local sources. The observation of 
the anisotropy of arrival directions of showers llT6l . llT9l 
supports this statement. At last, no dramatic fall of the 
flux of the primary particles relatively to the reference 
spectrum is seen at energies above 10^° eV. Indeed, due 
to lfT2l the number of expected events with energies 
above 6.3 • 10^^ eV is equal to 43.2 but only 13 events 
were observed. According to the reference spectrum the 
number of expected events is calculated as 16 that with 
the Poisson fluctuation taken into account agrees with 
the observed number. Due to ||T6| numbers of expected 
events with energies above 6 ■ 10^^ eV and lO^*' eV 
are equal to 167 and 35 accordingly while only 69 
and 1 events were observed. The reference spectrum 
gives 137 and 7 events accordingly. The disagreement 
decreased but not vanished. This fall may be regarded 
as the observation of the GZK supression of the flux of 
the primary particles. It may be commented that data 
lfT6^i include two empty bins which should contain 6 
events. The Poisson probability that no one was observed 
equals to ^ 2.5 • 10^'^. It should also be mentioned that 
the intensity of the reference spectrum is ^ 1.5 times 
higher than used in [161. Besides, some uncertainties 
in energy estimates are possible as we believe due to 
the constant integral intensity cut method [JJ which 
disregarded fluctuations. The data lfT3l illustrate possibly 
a contribution of local variable sources at energies above 
10^° eV while some uncertainties in energy estimates 
or aperture may be of importance at energies below 
5 • 10^^ eV. The same comments may be addressed to 
data IfTTl with additional remark that calculated energy 
estimate is 1.6 times less than used in data. Besides, 
data ifTSll . IfTTl show probably the second variable excess 
at energies (2-3)- 10^° eV. This second variable exess 
may be regarded as a contribution of heavy nuclei to 
the flux of the primary particles from the local sources. 
It is also not excluded that the Lorentz invariance may 
be violated at such huge energies lf20l . lf2Tl . If22l . The 
intensity of local sources on the Earth may be estimated 
as integral on differences J2 — Ji and J3 — Ji accordingly 
in case of the first exess. This integral is estimated as 
I K, A ■ 10^^^ TO~^s^^sr^^. If we assume a distance 
to local sources as i? ^ 30 Mpc, a typical energy 
E ^ 10^^ eV and an angle of emission as ~ 1 sr then 
we obtain the power of all local sources as 3 • 10'^^ W 
which is in agreement with estimates in jf23l . If24l . 

IV. Conclusion 

The new method has been suggested to estimate 
energy of extensive air showers by comparison all de- 
tector readings with calculated signals for a sample of 



individual events induced by various primary particles. 
Simulations of the individual shower development in 
the atmosphere have been carried out with the help 
of code CORSIKA-6.616 {?\ in terms of the models 
QGSJET2 |9| and Gheisha 2002 LIOJ with the weight 
parameter e = 10^^ (thinning). The program GEANT4 
1 11 1 has been used to estimate signals in the scintillation 
detectors from electrons, positrons, gammas and muons. 
New estimates of energy of the giant air shower observed 
at YA [71 have been calculated in terms of the QGSJET2 
1^ and Gheisha 2002 [101 models as £; - 2 • 10^" eV 
for the proton primaries and E ^ 1.7 ■ 10'^" eV for the 
primary iron nuclei. The base universal spectrum such 
as Jfc = A • £;-3-25 jj^yg been suggested at energies 
above lO" eV. It was also suggested that some possible 
local sources may produce variable contribution to the 
different regions of the energy spectrum at super high 
energies. The fall of the flux of the primary particles 
due to the GZK effect should be considered relatively 
the some base spectrum. Possibly, the second exess at 
energies (2-3)10^" eV has been observed in lEl, IfTTl . 
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